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Abstract 


Buckling is a Significant cause of structural failure. A 
theoretical investigation has been undertaken into the dynamic 
instability of complete spherical shells which are loaded 
impulsively and made from a viscoplastic material. . 


The critical mode numbers obtained in this thesis are 
Similar to those obtained for cylindrical shells having the same 
R/h ratios and material parameters while the threshold velocities 
were for higher viscoplastic spherical shells than for rigid- 
‘plastic spherical shells. 
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Introduction 


The analysis of dynamic buckling, of elastic-plastic 
and rigid plastic cylinders has been investigated by 
Abrahamson and Goodier [1] extended by Goodier [2], Goodier 
and Florence [3] and others [4-6]. Florence [7-8] and 
Wojewodzki [9-11] extended the analysis to viscoplastic 
cylindrical shells. Work on spherical shells has progressed 
through elastic-plastic and rigid plastic materials [12-14]. 
The analytical procedures of the above references were derived 
by assuming that buckling stemmed from the growth of small 
imperfections in otherwise uniform initial displacements 
and velocity fields. 

The dynamic buckling of a thin-walled complete spherical 
shell which is made from a viscoplastic material is 
investigated herein. A theoretical procedure is developed 
by first considering the unperturbed motion arising from 
a spherically symmetric initial velocity field. An 
examination is then made of the influence of small perturbations 
in this otherwise uniform velocity field and of the effect 
of small deivations from the initial sphericity of the shell. 
It transpires that certain harmonics grow very quickly and 
cause a spherical shell to exhibit a characteristic wrinkled 
Pee meeteeiseassumed that the Strains which are associated 
with the perturbed state are much smaller than those 
associated with the predominant motion. Thus, no unloading 


of the shell material occurs before the time when all the 





initial kinetic energy is absorbed by the spherical shell in 
the predominant mode. Moreover, it is further assumed that 
the wrinkled pattern in the shells is established prior to 


this time and is not modified by subsequent unloading. 
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Chapter 2 


Basic Equations 


The equilibrium equations, compatibility relations and 
strain and curvature relations have been presented in refs. [13] 
and [14] por a spherical shell. These results are repeated in 
this section for convenience. 

It may be shown for the particular case of a thin-walled 
spherical shell that the strain and curvature relations which 
were developed in reference [12] for an arbitrarily shaped 


shell, can be written 


¢ - 1 a 2 1 @ 
es = (Uy W + R We gg) (la) 
5 = {(cosecd Vv + cotd U - Ww oo it cosec’*¢ W W ) (1b) 
8 R cs) R as aos, 
2 = icosee U - cotd Vv + Vv “+ 
66 ° °#4R as "d 
2 cosecdo (W W + W W )] Gile) 
R ome ‘§g 'o 
en owe.) (1a) 
p R* “od 
2 = 1 2 2 2 2 
K 6 = aril + cosec’ ¢ Wr og + cotdo We 4) (le) 
: 2 : 
and K ye = pee occ Ce Wr og - cotd cosecd Wr 9] CIsE) 


when the displacements are defined as shown in Figure l. A 
straightforward modification of the general equations in 
ref. [12] for a shell of thickness h and density p gives the 


equilibrium equations 


= a ee 








Figure l: (a) Displacement and (b) Membrane Forces, Bending 


Moments and External Loads on a Spherical Shell. 
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N + cotd (Ny - No) + cosecd Nao - 


ohRU + P|R = 0 | (2a) 


o’ > 


cosec® Norg + Noorg * ZCOL® N50 - 
ohRV + PoR = 0 (2b) 
R 4 4 W, 4 W, 
and (Ny No) Ny (Ww 64 cotdé ») + 


2 
No cosec’@ Wr og + 
cosec¢ Nog (We 46 + Wr 94) ate 
(Nang We 4 + Wr oNgorg cosecd + 


Wo Noare coseco + 


cosec’$ Norg Wra) - (Mor og i 

2 cot Myig 7 cot? Morg + 

2M, + cosec*¢ Morag + 

2 cosecd me’ be + 

2 cosecd cotd Magra) - R2ohW + 

R°P, = 0 . (2c) 


2 
which are consistent with equations (la) and (1f) according to 
the principle of virtual work provided the influence of 
transverse shear deformation is disregarded. 
The strain and change of curvature relations (la) to (1f) 


Satisfy the compatibility equation 


= j= 





25 oh + cosec’¢ eee Eeoecot } =e = 
cots Bae a + aad - cosecd ©40'o0 = 
cosecd cotd © 40/9 + R(K 4 + Kg) = 0 (3) 


for infinitesimal displacements [12]. 

The constitutive equations for spherical shells which are 
Made from rigid viscoplastic material are developed in Appendix A. 
Equations (1.8a) to (1.8f) for a spherical shell made from rigid 


viscoplastic material can be written in the form 














kh : : 
Ny = = (2e, P Be) (4a) 
Y 
kh : : 
No = a (2e, =e e4) (4b) 
kh e 
Nog = — (4c) 
2 4 
kh ? . : 
M, = VU 2fe te pee (4d) 
12y* ¢ 2 
kh? : : 
Pee = (2k, + K,) (4e) 
7 12y* 8 b 
kh? «k 
ad Myg = ve | (A£) 
UZ eZ 
where k = Se) (5a) 
and y* = y(1 - k/VJ5) « (5b) 


Piewmetzniod Of Solution used here is Similar to that used in 
See eeeeez eancest ls} os) Tt 1s assumed that the displacement in a 


complete spherical shell consists of dominant parts 


= a> = 





(WwW, U = V = 0) and perturbation parts (w', U', V') so that 


W(t) + W'(o,8,T) (6a) 


W(¢o,0,T) = 
Wed, Oya) = UG, 0; TU) (6b) 
Bend V(6,8,T) = V'($,8,T). | (6c) 


Thus, the dominant strains and curvature changes are 


Prenat 
R 


e4 9 = 0. (FJa-c) 


se Kye = 0 (7d-£) 


and oe Se” 


a 


respectively, while the corresponding perturbation guantities 


are given by equation (la) to (lf) with primes. Similarly, 
Ny (0,9,T) = N (T) at Ny (0,9,7) (8a) 
Ng(¢,8,T) = N(t) + NA (6,8,T) (8b) 
Nog (0.9,T) = Nag 00 9-T) | (8c) 
My (O,9,T) = M(t) + My (,8,7) (8d) 
My ($,9,T) = M(t) + Ma (¢,8,T) (8e) 
and Mag (O.9,7) me sy ee (8£) 


If equations (6) and (8) are substituted into equation (2c), 
then the dominant and perturbation transverse equilibrium 


equations are 


(9a) 


N) 
ml al 
I 
N 
| 
ae 
ae) 
I 
= 
=) 
= 
© 


=: Mew 





and Rhee Ne) + N(W! + cotod W!, + cosec’®¢ W! = 


p “oo 


[Me oo + 2 €ocd Loo gs 


t 2 ' t 
2Mo + cosec’ ¢ Mar 96 + 2 cosecdo M4086 + 


$ 99) 


2 cosecd cotd Mier) ~ phR2W' = 0 (9b) 


respectively, when neglecting higher order quantities. 
It may be shown that in-plane equilibrium equations (2a) 


and (2b) are satisfied identically when 


1 


= —, [F + *o F + 
Ng 7 [ cosec’ é od cotd¢ Fy 4! (10a) 
i Si ee (lob) 
0 R2 od 
eel a 
and By) = ei cosecd [Fy 99 cotd F, 9] (1LGe) 
where F is a stress function and U = V = Eat = Po = 0. Now 


substituting the inverSions of equations (4a) to (4c) and 
eguations (l2a) to (12c) and (1d) to (le) into the compatibility 


equation (3) gives 


Bey 2) eet Rigger mest y? + 2)F = 0 (lla) 
U2 ae 2 _k _& 

where up [ khRe ' eo ae a mt (lib-e) 
and Weary =: cotd()r4 Heeocee O()),..- (11£) 


The transverse equilibrium equations (9a) and (9b) can be 


mer1tten in the form 


~u + 2(1 + a2)o0 + P = O (12a) 


oO TT 


= hr 





and Boge SSC es Se Sy ae 
ome 12(V- + 1) + Li¢(v2 + 2)u,, = 0 (12b) 


respectively, when using equations (1), (4), (6)-(8) and (10) 


and where 


and Bo =e = : (2c —c)) 


- 18 - 





Chapter 3 


Static Behavior 


Equation (12a) predicts 


: P 
= ae (13) 
2 
when a* << 1, while equations (lla) and (12b) become 
[2 (A, - 1) + 1] Or, - 2) Coe = 
3A, - a rs = 0 (14a) 
and oon mea oE Cine. ot: Oh Bilas 
eee) a Oe 2a 7 0 (14b) 
respectively, if it is assumed that 
u =a, P_(cos¢) 
Pi (cos$) (a. cos m@ + bon sin m0) (15a) 
and F = ConPry (cose) + 
Pi (cos$)(c_, cos m6 + d_ sin m8) (15b) 
where Ns = n(n + 1) (15c) 


and PL and Pe are Legendre and associated Legendre polynomials 
of degree n and order m. Two additional equations are also 
obtained but these are identical to equations (14a) and (14b) 


provide Dn and dan are substituted for aan ances. If ¢c 


n mn mn 


is eliminated from equations (14a) and (14b), then 


=< gee 











87wW a ot a = 0 (16a) 


~ wv a aaa ? = aa 

where vo, = + Oo (20, 3) O, Zee (16b) 
2r4_ - l 
n 
Now, substituting equation (13) into (16a) and integrating 
gives | E ng Cc | 
IR 

a =e + K (17) 
where K 1s a constant of integration which is related to the 
- initial imperfections in a shell. 8* which equals ae depends 
on strain rate as y* depends on strain rate as shown in 


appendix B. With a strain rate near zero, as would be the 


case for static bucklina 8* goes to infinity but 1/8* becomes 


=-5 ===. |= = ; ) (18) 


In the limit of En very small, ge << Y 





e 
limit L So 
e es = —~— € (19) 
a8 g 2 Cc e 
where Eo is the initial strain rate. Therefore equation (17) 
becomes ~ 4 
Sam eet] 
a =e n + K (20) 
mn 


when Eo IV 


This equation is of the form of creep buckling. 


= 2o- 








Chapter 4 


Dynamic Buckling 


4.1 Dominant Equation 


From equation (12a) 


Daag OG: (21) 


when a* << 1 and P= 0. oc is not time invariant and the 
solution of the above equation wn. ise done in two Sees, 

The first part will assume the value of ¢ to be constant with 
time and the second will give the exact solution, then the two 
solutions will be equated at the final displacements. For 


buckling the dominant stress is defined as 


i ho, (22) 


and equation (21) becomes 


e 
Voi te 273 —. (23) 
O 
O 
The solution of equation (23) is 
V T 
ui =—“27 fl- — (24a) 
C a4 
ie 
when satisfying the initial conditions W = 0 and W = We at 
t = 0 and where 
V O 
a oO mee (24b) 


is ‘the time when W ee 


a 





The exact solution uses equation (4a), (4b), (7a) and (7b) 


to revise o to reflect the time dependence as follows 


3 {- -!. (25) 


Me NE se Se CS 
= kh ~ 27 3 By Yor ee: 
where a ag a (27) 
e Reon 
Substituting equation (26), into (21) produces 
6C 
RY 


the equation for the exact solution. The solution of equation 
(28) is 
RY | 2V. 


———=—— 


= Meee: RY (29a) 


Ry 











e 


men Satisfies the initial conditions W = 0 and W = We at t = 


and 
RY ‘2V 
T. = — log, |1 + — 


£ (29b) 
6C RY 





is the time when W = 0. Setting U, (Te) for the first solution 
equal to Un (tT ,) for the exact solution and solving for Ce) Gn 
one obtains 

1 wiry)? | 2v 2V 


O 
—_ |— ——- log 1+ — 
Y3 IV Ry = RY (30) 


Q 


o lo 


De 





therefore, the dominant solution may be taken as equation (24a) 


with (24b) and (30) used with it. 
4.2 Perturbation Equation 
Using equations (15a-c) equation (12b) becomes 


N 
kh 


A. a + 


2 
BCA, = 2)e + wae 


ann’ tt mn 


2) AD = a a = 
= 6 [2 (A, 1) 1 OA, 2k one = 0 | 


Solving equation (14a) for Cron and substituting into 


equation (31) produces 








a = 
Qann’tt B O, Smn' st Ft wae Cnn ~ 0 
h 
where BS = n(n + 1) 
_ 3(A4_ - 2) 
rn 2 
er (2A, - 3) OY, - 2) 
4 te 
n 
Cc 1 Cc 
R? = = ao |. 
Ry* 2u Ry 


With a change of variable equation (32a) becomes 











‘a Cc Cc 
a; ~ TT, WwW + —/!ao_,-- 
Wbeletalca | aaa @| 2Ve Ry mm © 
: 3 3V 0 
T —+— ElrA a = 0 
28 5 Ry n mn 
where E =a — T/T ¢ 
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(31) 


(32a) 


(32b) 


(32c) 


(352d) 


(33a) 


(33b) 





iim 5) (33c) 


“G 
and u =) Se. (33d) 
C 


Putting equation (33a) into a format Similar to Florence [8], 


ene gets 








5 (2) 
a (1) “1H _ “7 af; 34 
Amn’ £E Qn ‘3 : mn’ é& ‘n mn (34a) 
C5 ae oes 
where oR T Ry OO (34b) 
(2) Cos 
QO. T Sais ome (34c) 
Oo 
3V 
and R? = <2 34 S ] ae =. (34d) 
n £ 2? Ry n 


As Re Maestrll a fUncCtiION Or time (&), in order to solve the 
equation uSing Florence's solution, the value of 3 + = | was 
set as a constant with respect to €. A value of 3/2 was used as 
we are primarily interested near € = 0 which negates the 


influence of the second term, thus 


mS 3/2 we (34e) 


z 
n n 
As given in Florence [8] the solution to equation (34a) 1s 


meee eae (2) GB (E)B. (35) 


<= Ae 





r(l + a. + b_) 
n n 


where Ann ‘© = = > a exp{- (ies, + 2) + ma) 
Nh 





ios 
(2) fe (Ge ne dL ae 38h) CC 


QUT (2) IM(Q, E) + t(AL +] + b,) M(2Q_) + 


2 M'(Q,) 5 H,0) (36a) 

r(i + a. + b,? 
and ee exp {—-(h. + 9.) + 
ae te H > 
n 
a3) 

Ae} (2) U(Q M(QE) - 
H(0,) 09,8) [26g (36b) 


In equations (36a) and (36b), M(Q_ &) and U (2, €) are shorthand 
notations to represent Kummer functions [15-17] (confuent 
hypergeometric functions) of the first and second kind, 


respectively, that is 


M(8_&) = M(l + a, + bo 2+ bos 2) (37a) 

U(2, &) = U(1l + a -- bie 2 + bo Q§). (375) 
Similarly 

M(Q,) =) 6 Bee an oF Du ose Doe Qe? (37¢) 

u(n,) EU(l+a +b, 2+b, 2) (374) 


= Ie a 





and primes over M(Q_) and U2) denote differential coefficients 








with respect to 2 evaluated at € = 1. The remaining symbols 
are 
9 (4) q(t) os 
ane 2 » boa (37e-£) 
2 Q 
n 
CE) 
2 9) 0 
oo fo) se ee || oe (37qg-h) 
n n n n ; 
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Chapter 5 


Peec mt omane pascuss ion 


Unfortunately, no experimental investigations appear to 
have been published on the dynamic instability of complete 
spherical shells and until reciently no theoretical work had 
been published. The results reported herein will then be 
compared with previous experimental investigations into the 
behavior of impulsively loaded cylindrical shells and in theory 
with the publications now in press. Equation (5.2) from 
Appendix E predicts mode numbers which are 3.173/3.624 = 0.876 
times as large as Jones and Ahn's [13] estimates for similar 
aluminum spherical shells and are eye = 1.08 times larger 
than Florence and Vaughan's [4] estimates for similar aluminum 
cylindrical shells. However, as shown in Figure 2, a plot of 
the first five sets of data from Table 1, the equation for the 
critical harmonic is 


605 


n = 2.0(R/h)° (38) 


which compares favorably with Lyons [6] curve with a slope of 


625. 
R/h Critical Harmonic n 
10 8 
bS..25 12 
30 a 
50 21 
100 32 


Table 1: A Comparison of the Critical Mode Numbers for a Visco- 
plastic Spherical Shell with Differing R/h values, = 
with R = 1.460 in., Vo = 9357.6 in/sec, y = 6500 sec 
oO, = 46,000 psi and p = 0.0002523 lb-sec*/in.”. 


a( D7 








10 20 30 50 100 


R/h 


Figure 2: Critical Mode of Buckling as a Function 
of R/h with we = 9357.58 in/sec and 


y = 6500 sec +. 


= Jen 





Figure 3 shows the growth of the amplitude function with regard 
to time for the continuous spectrum of the mode numbers n. 
Showing addition figures latér; this aqrowth is faster for 
larger ae larger y and for later staces of the process 

(small €). The change of preferred buckling mode is snoen 

by a dashed line. It is evident from Table 2a that equation (38) 
provides a reasonable estimate of the preferred mode of the 
amplification function which is obtained numerically from 

the corresponding theoretical analysis using Kummer's functions 
and numerical integration of the equations. The use of the 
Kummer's functions for solution ran into problems because of 
the size of the functions and that differences of these very 
large numbers were removing the first five or six significant 
figures. Even double precision was not eliminating the 
_-problem, for the asymptotic expansion [16] is a truncating 
series and even with double precision can't give the necessary 
Significant figures. 

The Florence method gave good results that correlated 
well, up to four significant figures, with the numerical 
method. The problem with the Florence method occurred when 
the values of inputs to the Kummer function became so large 
l.e. greater than 8 that one has to use an asymptotic 
expansion. Also the last couvle of computations with values 
Just less than 8 one loses some accuracy. As was observed 
by Florence [8] and Wojewodzki [11] the critical modes of 
viscoplastic cylinders were constantly lower than modes of 


Yrigid-plastic cylinders. Likewise as shown from Table 2a, 


eon = 





Ga OOS) A pue ‘cz°st = U/ae “*UT OOP'T = 4 
UFTM [TTeUuS TeoOTAeyds oTRSeTdoosTA e AOF SuoTtTzZoOoexZAsedut ARTOOTSeA (q) pue 


RZUowWeoeTdstq (2) xox FJ. QO°T pue g° ‘9° ‘7° 47° = 1 Fe UOTROUNY UOTReOTFITTdwy -¢ sANndbty 
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cers al Rese: (Koi) (1b-sec?/in.") ar 
Steel 1015 eae 102 .0007324 40.4 
AL 2024 T3 [1] 56 ~0002523 . 6500 
AL 6061 To [1] 44 .0002523 6500 
AL 6061 T6 [5] 44 -0902525 6500 
AL 6061 ™6 [6] 46 .0002523 6500 
Steel 1015 ‘wa 30 .0007324 40.4 


Table 2b: Values of Oy P and y that Correspond to Materials 


injnetaople Za. 


=. a 





Material Te C(10?in/sec) 
No. EQ. (24h) EO. (11d) 
1 094 ING AME 
2 110 16.3 
3 -108 16.3 
4 2085 16.0 
5 eL58 15.4 
6 - 140 ANS ees: 
a 2 EB2 55 
8 - 146 Be res. 
S, Shoe 135 
10 . 166 1s |S 
uae nat IS 13.6 
LEP Pelcy gk 13.6 
18 nOC7 Loe 2 
14 ~105 ike 
15 ae a 13.4 
16 2 ESO 13.4 
17 -142 eas 
18 5S S23 
19 © - 164. Loa! 
20 ~105 1 Seas 





Table 2c: Values of Te and C that Correspond to Materials 


in Table 2a. (Continued) 


Sees ee 








Material Tt C(10°in/sec) 


No. — BO. (24b) EQ. (114) 
21 120 ley 
22 .128 ony 
ze . 200 aA cal 
24 sige 14.0 
25 .169 13.9 
26 mie 14.0 
27 mek Sey 
28 .156 13.9 
29 .198 14.1 
30 .118 | 8.5 
31 mike i g.4 
32 | .145 8.6 
33 082 8.3 
34 071 8 
35 £157 8.7 
36 £119 8.5 
37 £115 8.5 
38 .158 8.7 
39 .103 8.4 
49 .087 8.3 
41 .076 8.2 


— 


Table 2c: Values of , and C that Correspond to Materials 


in Table 2a. 


a 3G 





the critical modes of spherical shells are considerably lower 
than those observed experimentally on cylindrical shells with 
the same material parameters and R/h ratios which were tested 
by the authors of references [1, 5, 6, 7], with the further 
observation that correlation seems to be closer at lower R/h 
values. Of course, n, in the present case, refers to the 
critical or most amplified harmonic. in the displacement 
field of a complete spherical shell which is expressed in terms 
of Legendre and associated Legendre polynomials (equation (15a)), 
while in reference [4] on cylindrical shells, it is associated 
with trignometric functions. 
The concept of threshold or cirtical impulse has been 
used by numerous investigators to mark the smallest impulse 
that a structure can tolerate without excessive deformation. 
It is evident from Figure 4 which is a plot of Table 3 that 
the displacement amplification function is a highly non-linear 
function of the impulsive velocity Oe The critical velocity 
is very dependent on the value of y used as shown in Figure 4. 
It 1s evident from Figure 4 that the displacement 
functions at T = Tp experience larger amplifications ee do 
the velocity functions. As stated by Jones and Ahn [13], this 
Suggests that shape imperfections might exercise a much more 
Significant influence on the instability of spherical shells 
than variations in the initial velocity field. . 


———— 
The critical harmonics are always rounded off to the nearest 
integer. 
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a a a SA ty typ aemespeeenennaseispegerenp-—ineern 


V Y Critical Critical 
. -1 Mode as Mode Bn 
(10 in/sec) sec ie n 


5 40 il Lei 2 0.12 
Mee 40 11 1.47 2 0.18 
ae 40 11 1.85 5 0.24 
11 40 11 2.05 5 0.26 
12 40 ia 2.28 6 ae 
15 40 11 Be7 7 0.40 
20 40 12 5.81 g 0.70 
30 40 13 paNG 10 2.38 
40 40 14 88.6 12 9.07 
A5 40 14 332. I NGG 
5 6500 Mi 1.60 2 0.19 
7S 6500 11 eS a 0.31 
10 6500 re 4.22 3 0.50 
11 6500 12 5.28 g 0.61 
12 6500 12 6.64 9 aS 
15 6500 12 me nG 10 eas 
20 6500 13 48.8 11 4.81 
25 6500 14 190. 12 17.4 
30 6500 15 788. ie 68.4 
B25 6500 15 1644. 14 138. 


Table 3: Critical Mode Numbers and Displacement Function 
7 Amplitudes Versus Initial Velocities at Differing 
Values of y with R = 1.460 in., R/h = 18.25, 


So = 46,000 psi and p = 0.0002523 lb-sec”/in.*. (Cont.) 


220m — 








V Y Critical | Critical 





(10 wees Cae Mode Avan Mode Bon 
n . : a n 
5 1,288,000 - Il 1.99 5 . Cass 
i) 1,288,000 us 3.85 8 0.48 
10 1,288,000 A 0 8.26 9 0.96 
11 1,288,000 12, 11.4 10 20 
12 1,288,000 13 Lo.9 10 1-76 
LS 1,288,000 a lies 45.7 11 4.73 
7.5 1,288,000 14 LV. 2 11.4 
20 1,288,000 14 304. 13 28.6 
25 1,288,009 ie) 2334. 14 205; 
5 900,000,000 iat 2.24 6 OS, 
Tees) 900,000,000 12 4.89 8 0.61 
10 900,000,000 12 TEU 10 1535 
11 900,000,000 3 i I Bs. 10 1,98 
12 900,000,000 13 257 iLat 2c) 
15 900,000,000 14 89.4 12 OR 9 
ere 900,000,000 14 274. 3 26.6 
20 900,000,000 15 856. 13 80.3 
22.5 900,000,000 LS Z5:02; 14 252% 
Table 3: Critical Mode Numbers and Displacement Function 


Amplitudes Versus Initial Velocities at Differing 
Values of y with R = 1.460 in., R/h = 18.25, 


o, = 46,000 psi and = 0.0002523 lb-sec’/in.*. 


- 40 - 





It should be noted that the present theory is an 
approximate one since the linearized Ageia of the constitutive 
equations for rate-sensitive material were used. In addition 
to the various simplifications and approximations which were 
introduced in the theoretical analysis, the influence of 
transverse shear stresses and deformations have been 
disregarded. Also, the influence of unloading has been 
neglected so that the effect of possible plastic behavior 
after T = Te has not been considered. Moreover, the influence 
of wave disturbance has been neglected. However, with final 
times (te) being in the neighborhood of 10-30 usec, which can 
be calculated by using equation (lle) and Table 2c, then one 
should not neaqlect wave phenomena. As an example materials 
1, 10, 20, 30, 40 have values of te = Seep bt 3 29.1, 


and 15.9 usec, respectively. 


ee ec 





Chapter 6 


Conclusions 


A theoretical investigation has been undertaken into the 
dynamic instability of complete spherical shells which are made 
from a viscoplastic material and loaded impulsively. 
Unfortunately, no experimental results are available to examine 
the validity, or otherwise, of the various theoretical predictions. 
However, the forms of the expressions for the critical mode 
number and threshold velocity are similar to those developed 
for cylindrical shells and spherical shells, except that the 
magnitude of the numerical coefficients are different. The 
threshold velocities for viscoplastic spherical shells came out 
higher than that for rigid plastic spherical seilile with the 
same R/h ratios and material parameters. The viscoplastic theory 
predicts critical mode numbers for spherical shells which are 
Similar to those observed experimentally on cylindrical shells 
with 10 < R/h < 30 but the correlation was better at the lower 
range the critical mode numbers were only one to three less 
than the rigid-plastic prediction where as the higher range 
produced mode numbers of four to eight less than the rigid 


plastic prediction. 





Chapter 7 


Recommendations 


To further the investigation of viscoplastic effects, the 
next logical step would be to look into the non-linear 
constitutive equations. This may result in a sufficient enough 


change so as to correct some of the problems of the threshold 


velocities. 
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Appendix A 


Constitutive Equations of Viscoplastic Spherical Shells 


The constitutive equations for cylindrical shells which 

are made from a rigid viscoplastic material are developed in 
e ® ¢ . 

references [8], [10], [11], and [18]. This appendix contains 
a brief derivation of constitutive equations for viscoplastic 
spherical shells. 

Consider a special linearized case of the constitutive 
equations for a rate-sensitive plastic material as derived by 


Perzyna [18] 
; oF 

ere. = 2kKY F 

tJ 00 


Cia) 





iJ 
where re and Gee are components of the plastic strain rate 
and stress tensors, k is the yield stress in simple shear and 
Y is a physical constant of the material later termed the 
viscosity constant. y is obtained by uniaxial tests as described 
later. In the above equation the function F, a function of 
strain rate and stress, is defined by F = V55/k - 1 where 
Jy = 1/2 od 4. Peto eso sever tonteome Lhe Stress deviator 
with components Si4° In accordance with the Levy-Mises flow 
rule for rigid perfectly plastic material, as y +®, J. = les 
where k = Aer 

The results of uniaxial stress experiments produced the 


Cowper and Symonds [19] power-type relation 


om 1/ 
Qs = Be 1+ s | Cl) 


mGAiG a= 





which give the non-linearized form of function F. Inverting 


(1.2) we get 


; O, 6) 
ge = yl—-il (1. 3a) 
e 
O 
Oo 
where O_ and —€ are the equivalent stress Tm eo. S ak and 
e e ; itis Re Big 
strain rate ie ees ae resvectively and o_ is the static 
3 13 #13 O 


uniaxial yield stress of the material and y and 6 are constants 


obtained experimentally for a given material 


o = [o Z = 6 O -2 ©) z + 30 24% (i331) 
e ible 11 ~22 22 eZ ° 

eee eae fF ca ae + c4,- £845 1° (eee) 
e jz 11 11 “$2 2D os ie 


Taking the linearized form of equation (1.2) we get 


E 
ce omen (1.4) 
e O Y 


where y is redefine as explained in Appendix RB. 


In the case of plane stress (053 = 053 = 032 = 0) 


equation (1.1) yields the following equations. 


> en ID One 922 
611.73 a I oe Ch Sa) 
VJ Zo - Oo 
boy = Lo 22 11 (SE) 
3 k Tp 
| Vaz ? 
ap tre = Dy 2 eae (1.5¢) 





Z 


where J. = i Oo 
ae 


, , k= oO /Y3. (1.5d-e) 


Inverting equations (1.5a-c) produces. 


O41 = ame 1d + E55) (1. 6a) 
ja (ee 166 
yo ~™ 22 il (1. ) 
See ae 
k dig 21. 
8] ee (1.6c) 
where y* = y(1 - kK/VJ5) = y(l - 0/5): (1.6d) 


Now, assume that the strain rates can be represented as follows: 


E14 =e) + Z Ky (1.7a) 
Eno = e5 + Z Ko (le 7b) 
E19 + E54 = C19 + Z Kio: (1.7c) 


Considering that y* is not constant and in fact is a function of 
Z, one assumes that the variation of y* with z is second order, 


a perturbation quantity. As shown in appendix B y* is dependent 


Ad 


° — 


on Ee, and as shown in Figure 5, e, = € 


o 1 ae : 
1 + a where Ey isa 
perturbation guantity dependent on z. Therefore, y* can be 
taken as a function of on anda mot a fUNCEITON Of Z. Now the 


constitutive equations can be integrated out as 
h/2 


(1.8a) 


Siow 





ane 
iy * 





olay 
Loe 





Dekh 
12y* 


= (400 














Eagubeso:y Variation of Strain Rate Across 


PhiGkhnheSs Om Section: 





Appendix B 


Variations of y and - 


Equation (1.4) from Appendix A is shown in figure 6 along 
with a straight line tangent to the Cowper and Symonds 
power type relation as suggested by Florence [7]. The equation 
Florence suggested 1S 


Ee = = Onl 


+ YE, (27215) 
where wv is the slope of the tangent line, tangent at the 


.* 
maximum initial strain rate Ear and o Pewee eIncencept OF 


O01 
this tangent line at static test strain rate. The slope of 
the tangent line 1s obtained by differentiations of equation (1.2) 


: -* 
and evaluating the results at E 





do Gis 

e Be Afsiee O 
Bee to eta /G (2.2) 

s * * 

de. i oy ex 


The value of o 1s now obtained by equation equation (2.1) and 


Ol 
mie 2eevaluatea at oa 


fe} = ©) 1+ j— l1-- (2.23) 


on" Oo 1 + = (22.45 





ese = 





where now To) and Oy ,/H compare with 6 and y from equation (1.4), 


respectively. Letting Oy )/H ae cy 


Yl 7s dk 
¥, = S€, ce 25 ar (2 35)) 
a 
Ee 
and oe To. 1+ —]. (2.6) 
Yq 


A replacement of - and y with So] and a respectively in the 
equations would give the identical solution with the result 


being the use of the tangent line approximation. 


The values of the material constants are 6 = 5 and 
, = 40.4 sa for hot rolled mild steel [20] and 6 = 4 and 
yY = 6500 sec * for 6061 T6 aluminum PlieeeenOm ein Sphnerica 


shells of radius R’the initial strain rate is computed as 


aa == 2V/R, where Ve is the uniform initial velocity. 


From Chapter 2 equation (5b) 


* 


Y = y(1 oe 06/9.) (2.7a) 


where k = a /¥3 and VI. = 0/3. (2.7b-c) 


Solving equation (1.4) for J 6/5 6 gives 


0 Y 
2 ene (2.8) 
Wes a sae § 


Since 1/y* is used more predominantly, substituting (2.8) into 
(2.7a) and inverting produces 


ee (2.9) 


Y Y Ee 


=) 





Figure 6: 





Cowper and Symonds Power Law and Linear 


Approximation to it [11]. 
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Appendix C 


Computer Proaqram Using Florence Method 


The following program was used to solve for mS) and 
Ban te): equations (36a) and (36b), by use of Mathematics 
Library Subroutine called HYPER which is a subroutine to solve 
Kummer's function of the first kind. Problems were encountered 
when taking the difference between functions as this was 


— 


eliminating the first five or six significant figures. 
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Appendix D 


Computer Program Using Numerical Method 


The following program solved the basic equation (34a) 
for pe) and Ban‘) by use of Mathematics Library Subroutine 
called DOSYK which is a subroutine to solve up to fourth order 
simultaneous differential equations with initial conditions. 


This program used equation (34d). 
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Appendix E 


Critical Impulse 


It 1S straightforward to show when — >> 1 that 1 can 
be eliminated from equations (34c) and (34e) and using 


equation (24b) to give 








eV. 2 8 R, $ 
ho << ees (Sa) 
Si /3 como) 
where 0/9. is defined in eguation (30). As shown in Table 4. 


the values of 0/5, lie within the range 0.98 < o/F, <b G 


(2) 


n 


with a medium value Oe ~ 1.07. Also the values of O lie 


(2) 


19 


within the range 1.27 < 0 < 4.42 with a medium value of 2.84. 
It is straightforward to show when using equation (34c) with 
n >> 1, that the critical harmonic 

n = 3.173(R/h) ? (5.2) 
(2) 


is associated with these values of oon and Q, : 
The threshold of dynamic plastic instability could be said 
to occur when Aan {9) = X, where X 1S an arbitrary positive 
number. If X = 100, then the initial displacement imperfections 
in the critical mode are magnified by a factor of one hundred. 


This amplification is achieved when Rie Pgeeaccording to 


Florence [3] asymptotic formula for equation (36a) 


iG Oe 





(A +1+b_)M(Q_)+2_M' (Q_) 
Aan (0) = exp |, + 2) | §———— (Ss 
oe Jake 


— 


where the right hand guantities are defined in equations 


i 7a-h) - 
(102insec) an Fo/% ae ae 
5 40 sO at 27, 
30 40 0.98 13 DDS 
5 6500 1.06 11 1.35 
30 6500 1.03 15 3.99 
5 1,288,000 1.11 11 1.41 
25 1,288,000 1.10 15 4.26 
5 900,000,000 1.16 il 1.47 
22.5 900,000,000 1.14 15 4.42 


(2) 


Table 4: Variation of 0/5, and Or 
and yY Wiehe — 1. 460ubhin kh, he — FE e725 oR = 46,000 psi 
and p = 0.0002523 lh-sec’*/in.”*. 


EOuvarclous ValLucs Of 


Finally equation (5.1) predicts a non-dimensionalized critical 


or threshold velocity 


="liil) h7R = (5.4) 





———— 
This asymptotic formula gave agreement with numerical values 
at —€ = 0 up to 4 significant figures. 


=G5a= 





This value does not compare favorably with Jones' [13] in which 
his coefficient is 17.8, since the velocities required should 


be higher than in the in rigid-plastic case. 
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